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Abstract. In this paper we study some optimization problems for nonlinear 
elastic membranes. More precisely, we consider the problem of optimizing the 
cost functional Sf(u) = Jqq f(x)u d'H N ~ 1 over some admissible class of loads 
/ where u is the (unique) solution to the problem — A p u + \u\ p ~ 2 u = in (1 
with \Vu\p~ 2 u v = / on dCl. 



In this paper we analyze the following optimization problem: Consider a smooth 
bounded domain ft C S. N and some class of admissible loads A. Then we want to 
maximize the cost functional 



Jan 

for / 6 A, where 7i d denotes the d— dimensional Hausdorff measure and u is the 
(unique) solution to the nonlinear membrane problem with load / 



Here, A p u = div(|Vu| p 2 Vu) is the usual p— Laplacian and ^ is the outer unit 
normal derivative. 

These types of optimization problems have been considered in the literature due 
to many applications in science and engineering, specially in the linear case p = 2. 
See for instance [5]. 

In recent years, models involving the p— Laplacian operator with nonlinear 
boundary conditions have been used in the theory of quasiregular and quasiconfor- 
mal mappings in Riemannian manifolds with boundary (see [9, 19]), non-Newtonian 
fluids, reaction diffusion problems, flow through porus media, nonlinear elasticity, 
glaciology, etc. (see [1, 2, 3, 8]). 

We want to stress that our results are new, even in the linear case. But since 
our arguments are mainly variational, and for the sake of completeness, we decided 
to present the paper in this generality. 

In this work, we have chosen three different classes of admissible functions A to 
work with. 

• The class of rearrangements of a given function fa. 

• The (unit) ball in some L q . 

• The class of characteristic functions of sets of given surface measure. 



1. Introduction 




N-l 



(1.1) 
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This latter case is what we believe is the most interesting one and where our 
main results are obtained. 

For each of these classes, we prove existence of a maximizing load (in the respec- 
tive class) and analyze properties of these maximizers. 

The approach to the class of rearrangements follows the lines of [6], where a 
similar problem was analyzed, namely, the maximization of the functional 



where u is the solution to —A p u = g in with Dirichlet boundary conditions. 

When we work in the unit ball of L q the problem becomes trivial and we explicitly 
find the (unique) maximizer for J, namely, the first eigenfunction of a Steklov-like 
nonlinear eigenvalue problem (see Section 4). 

Finally we arrive at the main part of the paper, namely, the class of characteristic 
functions of sets of given boundary measure. In order to work within this class, 
we first relax the problem and work with the weak* closure of the characteristic 
functions (i.e. bounded functions of given L 1 norm), prove existence of a maximizer 
within this relaxed class and then prove that this optimizer is in fact a characteristic 
function. Then, in order to analyze properties of this maximizer, we compute the 
first variation (or shape derivative) with respect to perturbations on the set where 
the characteristic function is supported. 

This approach for optimization problems has been used several times in the 
literature. Just to cite a few, see [7, 12, 15] and references therein. Also, our 
approach to the computation of the first variation borrows ideas from [13]. 

The paper is organized as follows. In Section 2 we include some preliminary 
results, some of which are well known but we choose to include them in order 
to make the paper self contained. In Section 3 we study the problem when the 
admissible class of loads A is the class of rearrangements of a given function /o. 
In Section 4, we study the simpler case when A is the unit ball in L q . Finally, in 
Section 5, we analyze the case where A is the class of characteristic functions of 
sets with given surface measure. 



In this section we collect some well known results that will be used throughout 
the paper. 

2.1. Results on rearrangements. First, we recall some well known facts on re- 
arrangements that will be needed in Section 3. 

Definition 2.1. Suppose f : (X, £,/i) — ► R+ and g : (X' , fj,') — ► K+ are mea- 
surable functions. We say f and g are rearrangements of each other if and only 



H({x G X: f(x) > a}) = fi'({x G X 1 : g(x) > a}), Va > 0. 

Now, given / G L'P(A), where A C R N with H d (A) < oo, the set of all 
rearrangements of f is denoted by TZf . Thus, for any / G TZf , we have 

H d {{x G A : f(x) > a}) = H d ({x G A : f Q (x) > a}), Va > 0. 
We will need the following Lemma, the proof of which can be found in [4] . 




2. Preliminaries 



if 
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Lemma 2.2. Let /o £ L p (d£l) and v e L p (dCt) such that fo,v > 0. Then there 
exists f G lZf such that 



[ fvdH N ~ 1 = sup f hvdH N - : 
Jan hen Sa Jdn 



The following result can be easily deduced from [17] (Theorem 1.14 p. 28). 

Theorem 2.3 (Bathtub Principle). Let (f2,E,/i) be a measurable space and let f 
be a real-valued, measurable function on Q such that [i({x : f(x) > t}) is finite for 
all f £ I. Let the number G > be given and define the class C of measurable 
functions on f2 by 

C = {g: < g(x) < 1 for all x and / g(x) dfi — G}. 

Jn 

Then the maximization problem 



I = sup / f(x)g(x)d(i 
gee Jn 

is solved by 

(2-1) g(x) = X{f>s}(x) + c X {f= s }(x), 

where 

s = mi{t : n({f > t}) < G} 

and 

cv({f = s})=G-»({f>s}). 
The maximizer given in (2.1) is unique if G — fi({f > s}) or if G = m({/ > s })- 

2.2. Results on differential geometry. Now we state without proof some results 
on differential geometry that will be used in the last section. The proof of these 
results can be found, for instance, in [14]. 

Definition 2.4 (Definition of the tangential Jacobian). Let £1 C M. N be a smooth 
open set ofR N . Let $ be a C 1 field overM. N . We call the tangential Jacobian of $ 

J T ($) := H*']" 1 "!^*). 

where v is the outer unit normal vector to dfi, <&' denotes the differential matrix of 
J(3>) is the usual Jacobian of $ and T A is the transpose of the matrix A. 

The definition of the tangential Jacobian is suited to state the following change 
of variables formula 

Proposition 2.5. Let f e L 1 {^{dVL)). Then f o $ e L 1 (dQ) and 



f fdH N -'= f (fo^)J T ^)dH N ' 

J<S>(dn) JdQ 



Definition 2.6 (Definition of the tangential divergence). Let W be a C 1 vector 
field defined on M. N . The tangential divergence of W over dfl is defined as 

div T W := divW- (WV, v), 

where v is the outer unit normal vector to dfl and (•, •) is the usual scalar product 
in R N . 
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With these definitions, we have the following version of the divergence Theorem. 

Theorem 2.7. Let SI be a bounded smooth open set ofR N , D C <9S1 be a (relatively) 
open smooth set. Let W be a [W 1 ' 1 (dil)] N vector field. Then 

f divrWdW 1 * -1 = / (W,^ T )dH N - 2 + f H(W,i/)dW N -\ 
Jd Jon Jd 

where v T is the outer unit normal vector to D along <9S1 and H is the mean curvature 

of on. 

3. Maximizing in the class of rearrangements 

Given a domain SI C M. N (bounded, connected, with smooth boundary), first we 
want to study the following problem 



(3.1) 



j-A p u+ \u\p- 2 u = in SI, 

|| Vu |P-2|m = / Qndn . 

Here p e (1, oo), A p u = div(|Vu| p_2 Vu) is the usual p— Laplacian, is the outer 
normal derivative and / G L q (dil) with q > jp . 
We say u <G M /1 ' P (S1) is a weak solution of (3.1) if 

/ \Vu\ p - 2 \7uVv + \u\ p - 2 uvdH N = [ fvdH 1 *- 1 
Jn Jon 

for all v e W 1 *^). 

The restriction q > jp- is related to the fact that = where = p(iV — 
1)/(7V — p) is the critical exponent in the Sobolev trace imbedding W 1,P {Q) 
L r (dil). So, in order for that the right side of last equality to make sense for 
/ e L q {dVL) we need w to belong to L q (SI). This is achieved by the restriction 
q' < P*. 

It is a standard result that (3.1) has a unique weak solution Uf, for which the 
following equations hold 

(3.2) f fu f m N ~ 1 = sup l(u), 

Jdn uew 1 .p{n) 

where 

l(u) = -^—{p ( fudH*- 1 - { \Vu\ p + \u\ p dH N ). 
V ~ 1 L Jdn Jn > 

Let fo S L 9 (<9S1), with q = p/(p— 1), and let lZf be the class of rearrangements 
of fo. We are interested in finding 



(3.3) sup / fufdH 1 *- 1 . 

feiz fo Jon 

Theorem 3.1. There exists f £ lZf such that 

J(f)= f fudU N - 1 = sup J(f)= sup f fu f dH N -\ 
Jon f£K fn fen fn Jon 



where u — Uf. 
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Proof. Let 

/= sup f fufdH 1 *- 1 . 
fen Sa Jan 

We first show that 7 is finite. Let / G 7£/ . By Holder's inequality and the trace 
embedding we have 

/ |Vu / | 1 ' + | U/ | 1 'dW JV <C||/|| i , { 8 n )||« / || w i.p( n ), 
Jn 

then 

(3.4) \\u f \\ W i, P{ n)<c yfen fo 

since \\f\\ L i(dn) = ||/o ||x,«(an) for all / e TZ fo . Therefore 7 is finite. 

Now, let {fi}i>i be a maximizing sequence and let m — Uf t . From (3.4) it is 
clear that {iti}i>i is bounded in W 1,P {Q), then there exists a function u e W 1,P (Q) 
such that, for a subsequence that we still call {ui}, 

Ui — 1 u weakly in W 1,P (VL), 
Ui — > u strongly in L p (fl), 
Ui — > u strongly in L p (dfl). 

On the other hand, since {fi}i>i is bounded in L q (dQ), we may choose a subse- 
quence, still denoted by {fi}i>i, and / G L q (dSl) such that 

fi f weakly in L 9 (dQ). 



Then 



7 = lim / fiUidH 
'an 



i — *oo 



JV-1 



= J— Km \ p f f iUidn "-i- f \Vu t \ p + \ Ul \ p dH N \ 

< ~^—\pf fu&H"- 1 - ( \Vu\ p + \u\ p dU N \. 
P - 1 L Jan ' Jn - 1 

Furthermore, by Lemma 2.2, there exists / £ 7£/ such that 

/ /udH^- 1 < / /udW 7 
Jan Jan 



/JV-l 



Thus 



7 < -J— fp / fudU 1 *- 1 - [ \Vu\ p + \u\ p dH N \. 
P - 1 L Jan Jn > 

As a consequence of (3.2), we have that 

7 < -^—{p\ fudH 1 "- 1 - { \Vu\ p + \u\ p dH N ) 
P - 1 Wan 7o J 

< -^—{p I fudU 1 "- 1 - I' \Vu\ p + \u\ p dU N ) 
P - 1 Wan Jn J 



AT-l 



= / fudH 

Jan 
< 7. 

Recall that u = Uf. Therefore / is a solution to (3.3). This completes the proof. □ 
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Remark 3.2. With a similar proof we can prove a slighter stronger result. Namely, 
we can consider the functional 

Mf,g)-= [ gudH N + [ fudH N -\ 

Jn Jon 
where u is the (unique, weak) solution to 

{—A p u + \u\ p ~ 2 u = g in 0, 
W~ 2 ^=f on 

and consider the problem of maximizing J\ over the class lZ go xTZj„ for some fixed 
.90 and f . 

We leave the details to the reader. 



4. Maximizing in the unit ball of L q 
In this section we consider the optimization problem 

maxj(f) 

where the maximum is taken over the unit ball in L q (dil). 

In this case, the answer is simple and we find that the maximizer can be computed 
explicitly in terms of the extremal of the Sobolev trace embedding. 

So, wc let / £ L q (d£l), with q > jp, and \\.f\\Li(dn) < 1) we consider the problem 

(4.1) sup f fu s m N , 

where Uf is the weak solution of 

f-A„u+ \u\p- 2 u = in Q, 



(4.2) 



.pi. 

| Vu |p-2&u = f Qn g n 



dv 

The restriction q > jp is the same as in the previous section. 

In this case it is easy to see that the solution becomes / = v^ 1 where v q i £ 
W 1,P (Q) is a nonnegative extremal for S q i normalized such that ||?y llz^'ran) = 1 
and S q i is the Sobolev trace constant given by 

Si= . nf J n \Vv\ p + \v\ p dH N 
q ' vew-nn) (j gn \ v \ q > dn N-i}f 

Furthermore u = Ui = 1 h_ 1 v q i. Observe that, as q' < p* there exists an extremal 

Si 

for S q i. See [11] and references therein. 
In fact 



J(f)= [ fu<M N - 1 = f \Vu\ p + \u\ p dH N 
Jan Jn 



\Vv„i\ p + \v n i\ p dH 



N 



«P/(P-1) J n l q n 1 q 1 " Cl/(P-1)' 

q' " q' 
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On the other hand, given / £ L q (dQ), such that < 1, we have 

J(f)= [ /u/dW^- 1 < ||/IU.(an)||«/||^ ( an) 



g 

from which it follows that 

This completes the characterization of the optimal load in this case. 

5. Maximizing in L°° 
Now we consider the problem 

N-l 



(5.1) sup / 4>Uj,d?i 

0GB Jan 

where B := {<j> : < <p(x) < 1 for all x E dQ and j gn cjxM 1 *- 1 = A}, for some 
fixed < A < H N ^ 1 (dQ), and is the weak solution of 



(5.2) 



-A p m+ \u\ p - 2 u = infi, 



This is the most interesting case considered in this paper. 

5.1. Existence of optimal configurations. In this case, we have the following 
theorem: 

Theorem 5.1. There exists D C <9f2 with H N ^ 1 (D) = A such that 

f XDUodH"- 1 = sup / (fm^dH N -\ 
Jon 0eB Jan 

where Ud = u XD . 
Proof. Let 

/= sup f <jm^m N - 1 . 

Arguing as in the first part of the proof for Theorem 3.1 we have that I is finite. 

Next, let {4>i\i>\ be a maximizing sequence and let u, = u^ >i . It is clear that 
{ui}i>i is bounded in W 1,p (fl), then there exists a function u £ VF 1,p (f2) such that, 
for a subsequence that we still call {ui\i>\ 

u l — u weakly in W 1,P (Q,), 
Ui — > w strongly in L p (f2), 
Uj — > it strongly in L p (dfl). 

On the other hand, since {0i}i>i is bounded in L°°(dCl), we may choose a subse- 
quence, again denoted {</>i}i>i, and £ L°°{d£l) and such that 

c^ A 4> weakly* in L°°(dSl). 
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Then 



lim / 



an 



Ui dH N ~ l 



= -J—\i m \pf famm 1 *- 1 - [ \v Ui \ p + \ Ui \ p dH N } 

V - 1 L Jan Jn > 

< ~^—\pf cj>uAU N - 1 - { \Vu\ p + \u\ p dH N ). 
P - 1 L Jan Jn > 

Furthermore, by Theorem 2.3, there exists D C dVL with H N ^ 1 (D) — A such that 



/ (^udH*- 1 < I XDudH N -\ 
Jan Jan 

and 

{t < u} C D C {t < u}, t := inf{s : H N -\{s < u}) < A}. 

Thus 

/ < -^—{p I XDudH"- 1 - [ \X7uf + \u\PdH N \. 
P ~ 1 L Jan Jn > 

As a consequence of (3.2), we have that 

-^{p [ XDudH*- 1 - ( \Vu\ p + \u\v&H N X 
- 1 L Jan Jn ' 

^{pJ^XDUodH"- 1 - jf \Vu d \p + \u D \ p dH N } 



I < 



< 



las 



N-l 



Xdu d dH 

I an 
< I. 

Recall that ud = u XD . Therefore xd is a solution to (5.1). This completes the 
proof. □ 

Remark 5.2. Note that in arguments in the proof of Theorem 5.1, using again the 
Theorem 2.3, we can prove that 

{t < u D } C D C {t < u D } 

where t := inf{s : 7Y Ar_1 ({s < ud}) < ^4}- Therefore ud is constant on 3D. 

5.2. Domain Derivative. In this subsection we compute the shape derivative of 
the functional J{xd) with respect to perturbations on the set D. We will consider 
regular perturbations and assume that the set D is a smooth subset of dfl. 

Then, by using the formula for the shape derivative, we deduce some necessary 
conditions on a (regular) set D in order for it to be optimal for J in the L°° setting. 

Also, this formula could be used to derive algorithms in order to compute the 
actual optimal set (cf. with [10]). 

For the computation of the shape derivative, we use some ideas from [13]. 

We begin by describing the kind of variations that we are considering on the set 
D. Let V be a regular (smooth) vector field, globally Lipschitz, with support in a 
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neighborhood of dfl such that (V, v) = and let ip t ■ K w — ► 1* be defined as the 
unique solution to 



V>o(z)=z xeR N . 



(5.3) 

We have 

i>t(x)=x + tV{x) + o(t) VxeR N . 
Now, if D C dCl, we define D t := ^ t (£>) C 90. 

First, we compute the derivative at t = of the surface measure of the set D t . 
That is, we want to compute 

dt t=o 

Lemma 5.3. With the previous notation, if D C dQ is a smooth (relatively) open 
set, then 



i*»-(A) 



= I dh 
-=o J D 



divVdH 



N-l 



Proof. We will use the following asymptotic formulae, for which the proofs can be 
found in [14]: 

(5.4) J^(*) = l + *divV(x) + o(t), 

(5.5) [ipt 1 ]'(x)=Id-tV(x) + o(t). 

Then we have, by the change of variable formula, Proposition 2.5, 

H N -\D t ) = / dW"-^ / \[^ 1 }'(x)iy\JMx)dn N - 1 . 

J D t JD 

Hence by (5.4), (5.5) and the definition of J T we get, using that (V, v) = 0, 
H N - 1 (D t )^n N - 1 (D)+t f divV dH 1 *- 1 +o(t). 

JD 

Therefore, we arrive at 
d 



f t H»->(D t ) 



divVdH*- 1 

t=o 



This is what we wanted to show. □ 
Now, let 



(5.6) 



I(t) = f u tX D t dH N -\ 
Jon 

where u t G W 1,P (Q) is the unique solution to 

f-ApUt + \u t \P- 2 u t = inQ, 
\\y Ut \p-^ = XDt ondQ 

and assume that D C dSl is again a smooth (relatively) open set. 
We have the following Lemma: 

Lemma 5.4. Let uq and u t be the solution of (5.6) with t = and t > 0, respec- 
tively. Then 

u t u in W 1,p (fi), as t 0+. 
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Proof. The proof follows exactly as the one in Lemma 4.2 in [6]. The only difference 
being that we use the trace inequality instead of the Poincare inequality. □ 

Remark 5.5. It is easy to see that, as ip t ~ > Id in the C 1 topology, then from 
Lemma 5.4 it follows that 

wt := Ut o ip t — > mo strongly in W 1,p (il). 

Now, we arrive at the main result of the section. 

Theorem 5.6. With the previous notation, if D C dil is a smooth (relatively) open 
set, we have that I(t) is differentiate at t = and 



±I(t) =—J u (V,v r )dH N -' 
dt t=o p-1 J dD 



where u is the solution of (5.6) with t = and v T stands for the exterior unit 
normal vector to D along dfl. 

Proof. By (3.2) we have that 

I(t)= sup -^—{p[ vxDtdH"- 1 - [ \Vv\P + \v\PdH N \. 
vew 1 ^^) P - 1 I Jan Jn J 

Given v G W 1,p (Cl) we consider u = v o tjj t e W 1,p (Cl), then, by the change of 
variables formula, Proposition 2.5, 

f vxDtd^-^f u X DJrAdH N - 1 
Jan Jan 

= [ ux D dH N - 1 +t[ uxodivrVdH 1 *- 1 + o(t). 
Jan Jan 

Also, by the usual change of variables formula, we have 

\Vv\ p dH N = [ \ T W t \- 1 {x)Vu T \ p J^ t dH N 
Jn 

= f \(I -t T V' + o(t))Vu T \ p {l + tdivV + o{t)}dH N 
Jn 

= [ {\Vu\ p -tp\Vu\ p - 2 (Vu, T V'Vu T )+o(t)}{l + tdivV + o{t)}dU 
Jn 



Vu\ p dH N + t / \Vu\ p dWVdH N 
n Jn 

tp f \Vu\ p - 2 {Vu, T V'Vu T )dU N + o(t), 
Jn 



and 



f \v\ p dH N = [ \u\ p J?p t dH N = [ \u\ p dH N +t [ \u\ p divVdH N + o(t). 
Jn Jn Jn Jn 
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Then, for all v G W 1 ' P (Q) we have that 

p [ vxdAK"- 1 - I \Vv\ p + \v\ p dH N 
Jon Jn 
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=p [ uxodH*- 1 - [ \Vu\ p + \u\ p dH 
Jdn Jn 

I 



'on 
+ t 



p I uxd&™ t VAU n - / {\Vu\ p + \u\ p )diWdH 
on Jn 



+ o(t). 



N 



+ P f \Vu\P- 2 (Vu, T V'Vu T )dH T 
Jn 

Therefore, we can rewrite I(t) as 

I{t) = sup -^—{ip{u) + t<t>{u) + o(t)}, 
uewi.p(n) p—l 



where 



and 



<p(u)=p f uxodH"- 1 - f \Vu\ p + \u\ p dH N 
Jdn Jn 



0(u) =p [ u XD div T V dH 1 *- 1 - I (\X7u\ p + \u\ p )divVdH N 
Jan 



+ P f \Vu\ p - 2 (Vu, T V'Vu T )dH N . 
Jn 



If we define Wt — u t o ip t for all t we have that wq = uq and 

1 



I(t) = —j{^( Wt ) + t<t>{w t ) + o(t)} 



for all t. Thus 



J(t) - 1(0) > -^{^(ii ) + t<j>{u ) + o(t)} - ^yvM, 



then 
(5.7) 

On the other hand 



liminf ^)-/(Q) > i M 
*^o+ t ~ p- r y °' 



J(t) - 7(0) < -J— {p(tu t ) + t0(«; t ) + o(t)} - -J— V («; t ), 
p — 1 p — 1 



hence, 



By Remark 5.5, 

therefore, 
(5.8) 



J(i) -7(0) 1 ,. , 1 . . 



<p(w t ) -> 0(t*o) as f -> + , 
J(t) - 7(0) _^ 1 



lim sup 

t-*o+ 



< —^(uo). 
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From (5.7) and (5.8) we deduced that there exists I'(0) and 

j'(0) = ^-j4>(u ) 

= -^{p f uoXD<Uv T VdH N - 1 +p f \Vu \ p - 2 (Vu Q , T V'VuZ)dH N 
P ~ 1 I Jan Jn 

-J (|Vu | p + |u | p )divV'dW Jv j. 

Now we try to find a more explicit formula for I'(0). 

In the course of the computations, we require the solution u to 

J -Au + \uo\ p ~ 2 u = in 0, 
[\Vu \p-^= Xd ondQ, 

to be C 2 . However, this is not true. As it is well known (see, for instance, [19]), u 
belongs to the class C 1 ' 6 for some < 6 < 1. 

In order to overcome this difficulty, we proceed as follows. We consider the 
regularized problems 



(5.9) 



'-div((|Vu§| 2 +£ 2 )<P- 2 )/ 2 V^) + \u%\p- 2 u% = in n, 
(|V U g| 2 + £ 2 )(f- 2 )/ 2 ^=XD on 30. 



It is well known that the solution ug to (5.9) is of class C 2,p for some < p < 1 
(see [16]). 

Then, we can perform all of our computations with the functions Uq and pass to 
the limit as e — > 0+ at the end. 

We have chosen to work formally with the function u in order to make our 
arguments more transparent and leave the details to the reader. For a similar 
approach, see [13]. 



Now, since 



div(K|no = p\u \p- 2 u (\7u ,v) + Kl p chw, 

div(|Vu | p V) = p\Vu \ p - 2 (Vu D 2 u ,V) + \Vu \ p dWV, 



we obtain 



I'(0) = -^—\pf UoXodWrVdH^+p f \Vu \ p - 2 (Vuo, T V'Vv%}dH 
P ~ 1 I Jan Jn 

f div((\X7u \ p + \u \ p )V)dH N +p [ |V U o| p " 2 (V Uo ^ 2 Mo,^}dH / 
Jn Jn 

+p f \u a \ p - 2 u (Vu ,V)dH N 
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Hence, using that (V, v) = in the right hand side of the above equality we find 

I'(0) = -M / u oXndiv T V dH"- 1 
P - 1 I J an 

+ [ \yuo\ p - 2 (Wu Q , T VVuS + D 2 u V T )dn N 
Jn 

+ jf luoMuo^Vuo.VOdW^j 
= -^r( / «oXDdiv T ydW JV - 1 + / |Vu r 2 (Vu ,V((Vuo,^)))dW JV 



n 



Since uq is a week solution of (5.6) with t — we have 

/'(o) = / u oX Ddw T vdn N - 1 + [ (Vu ,v)xDdn N - 1 \ 

P - 1 I Jan Jen J 

= / div T (u V) X Dd'H N - 1 

P - 1 Jan 



p-1 



P ' u (^,^)dW JV - 2 . 



ar> 



This completes the proof. □ 

The following corollary is a result that we have already observed, actually under 
weaker assumptions on D, in Remark 5.2. 

Nevertheless, we have chosen to include this remark as a direct application of 
the Lemma 5.3 and Theorem 5.6. 

Corollary 5.7. Let \d be a maximizer for J over the class B and assume that 
D C dSl is a smooth (relatively) open set. Let up be the solution to the associated 
state equation 



-A p u + \u\ p ~ 2 u = in SI 

<du 
dv 

Then, u D is constant along 3D. 



N-2 



Proof. Recalling the formula for the derivative of the volume, that is, 

±H N -\D t ) = f div T VdH N - 1 = [ (V,v T )dH 
dt t=o J D J dD 

and the fact that D is a critical point of /, we derive 

d 



i'(o) = c gn N - 1 (D t ) 



^=^> u = constant, on dD. 
t=o 



As we wanted to prove. □ 
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5.3. Final comments. It would be interesting to say more about optimal config- 
urations. For instance: 

• What is the topology of optimal sets? Are optimal sets connected? 

• What about the regularity of optimal sets? Is it true that the boundary of 
optimal sets are regular surfaces? 

• Where are the optimal sets located? 

These questions, we believe, are difficult ones and we can only give an answer 
in the trivial case where the domain f2 is a ball. In this case, by symmetrization 
arguments (by means of the spherical symmetrization, cf. with [12, 18]) it is straight 
forward to check that optimal sets are spherical caps. 

This example also shows that the uniqueness problem is far from obvious. 
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